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Abstract 

We provide complete analytic expressions for the isolated prompt photon 
production cross section in e + e~ annihilation reactions through one- loop or- 
der in quantum chromo dynamics (QCD) perturbation theory. Functional 
dependences on the isolation cone size 5 and isolation energy parameter e are 
derived. The energy dependence as well as the full angular dependence of the 
cross section on 6> 7 are displayed, where # 7 specifies the direction of the pho- 
ton with respect to the e + e~ collision axis. We point out that conventional 
perturbative QCD factorization breaks down for isolated photon production 
in e + e~ annihilation reactions in a specific region of phase space. We discuss 
the implications of this breakdown for the extraction of fragmentation func- 
tions from e + e~ annihilation data and for computations of prompt photon 
production in hadron-hadron reactions. 
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I. INTRODUCTION 



The cross section for the inclusive yield of high energy photons in hadronic final states 
of e + e~ annihilation is well-defined and, as demonstrated in our earlier paper [|l], it may 
be calculated reliably within the context of quantum chromodynamics (QCD) perturbation 
theory. However, an important practical limitation of high energy investigations, whether 
in hadron-hadron reactions or e + e~ annihilation processes, is that photons are observed and 
their cross sections are measured only when the photons are relatively isolated, separated 
to some extent in phase space from accompanying hadrons. Experimental procedures vary, 
but the essence of isolation is that a cone of half-angle S is drawn about the direction of 
the photon's momentum, as sketched in Fig. [I], and the cross section is defined for photons 
accompanied by less than a specified amount of hadronic energy in the cone, e.g., E^ ne < 
E max = e^-Ky. In this paper, we undertake a thorough analysis of isolated prompt photon 
production in e + e~ annihilation. Previous theoretical studies of isolated prompt photon 
production in e + e~ — > 7X include those of Refs. 

In this paper, our calculations of photon yields in e + e _ — ► 7X are carried out through 
one-loop order. We compute explicitly direct photon production through first order in 
the electromagnetic coupling strength, a em , and the quark-to-photon and gluon-to-photon 
fragmentation contributions through first order in the strong coupling strength a s . We 
display the full angular dependence of the cross sections, separated into longitudinal sin 2 9 1 
and transverse components (1 + cos 2 7 ), where # 7 is the direction of the 7 with respect to 
the e + e~ collision axis. 

A proper theoretical treatment of the isolated energetic photon yield requires careful 
consideration of the origins of both infrared and collinear singularities in QCD perturbation 
theory. In a theoretical calculation, photon isolation limits the final-state phase space ac- 
cessible to accompanying gluons (g) and quarks (q). This phase space restriction breaks the 
perfect cancellation of soft singularities in each order of perturbation theory that guarantees 
reliable predictions in the inclusive case ||. An uncanceled infrared singularity signals a 
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breakdown of conventional perturbative factorization. In this case, the uncanceled infrared 
singularity leads to an inverse-power divergence at the level of the partonic cross section 
[i.e., a 1/(1 — xi) divergence as x± — > 1] in the order a s one-loop quark-to-photon frag- 
mentation term. After convolution with the parton-to-photon fragmentation function, this 
inverse-power divergence becomes a logarithmic divergence in the isolated cross section. 

Owing to isolation, the predicted photon cross section develops explicit functional de- 
pendence on the isolation parameters and 5. We show that at one-loop level, this cross 
section is well defined in most of phase space, so long as and 5 are not too small. However, 
the region near x 1 = 1/(1 + e^) is special, and the perturbatively calculated cross section is 
singular there; x 7 = 2E 7 / -y/i, where a/s denotes the center-of-mass energy of the e + e~ reac- 
tion. In earlier papers |6j , we presented analytic expressions for the dependence on isolation 
parameters in the specific case of photon production in hadron-hadron collisions [0-0. In 
this paper, we provide a much more detailed treatment of photon isolation, concentrating 
on energetic photon production in electron-positron annihilation: e + e~ — > 7X @||. 

In a perturbative QCD (pQCD) calculation of the inclusive yield of photons, the quark- 
photon collinear singularities that arise in each order of perturbation theory, associated with 
the hadronic component of the photon, are subtracted and absorbed into quark-to-photon 
and gluon-to-photon fragmentation functions, D(z, fi 2 ), in accord with the factorization the- 
orem [ 1Q| . The scale fi 2 denotes the fragmentation scale that separates the non-perturbative 



domain from the region in which perturbation theory should apply; z is the momentum frac- 
tion carried by the observed photon from its parent parton. Since fragmentation is a process 
in which photons are part of quark, antiquark, or gluon "jets", it is evident that photon 
isolation reduces the contribution from fragmentation terms. For a small value of the energy 
resolution parameter, e^, the isolation cut eliminates most of the contribution from parton- 
to-photon fragmentation. However, may never be equal to zero either experimentally or 
theoretically. Experimentally, the finiteness of is guaranteed by detector resolution. The- 
oretically, the perturbative calculation of the cross section for isolated photons is ill defined 
if €h vanishes |§. 



Fragmentation is modeled in perturbation theory as a collinear process, whereas experi- 
mental cone sizes are finite (5 7^ 0) and partons are manifested as sprays of hadrons. Cor- 
respondingly, there is an inherent conceptual incompatibility between theoretical, collinear 
fragmentation functions and empirical fragmentation functions, D exp (z, /i 2 ,<5), that more 
naturally would be defined with reference to a cone of specified size. In this paper, we 
limit ourselves to the usual collinear fragmentation functions, D c ^ (z, /i 2 ). We derive the 
expected dependence of the cross section on cone size 5 and energy isolation e^, and we 
display a limited region of phase space in which the incompatibility of the collinear and fi- 
nite cone size assumptions leads to collinear sensitivity of the isolated photon cross section. 
We discuss the impact of both infrared and collinear sensitivity on computations of prompt 
photon production in electron-positron and hadron-hadron reactions. 

All four groups at LEP have published papers on prompt photon production [|TTJ. A 



measurement of the photon fragmentation function from an analysis of the two-jet rate in 
Z decays is reported by the ALEPH collaboration fll2| . Isolated prompt photon production 
has been treated theoretically previously, but our approach is different. In Ref. 0, the 
authors concentrate on events having a specific exclusive topology, such as a photon plus 
one hadronic jet; they discuss the extraction of the quark to photon fragmentation function 
from such data. In Ref. |§, isolated direct photon production through order a em is treated; 
the fragmentation terms there are included only at lowest order, not through order a s , and 
the angular dependence of the cross section is not derived. Practical aspects of confronting 



theoretical calculations with data from LEP are addressed in Ref. JT3[ . 

In Section |T[ we discuss the similarities and differences in the theoretical calculations 
of the cross sections for inclusive and isolated photons, and we provide a general factorized 
formula for the isolated photon cross section. Based on the definition of isolation, it is clear 
that the isolated photon cross section is part of the inclusive cross section. We write the 
isolated cross section as the inclusive cross section minus a subtraction term IHI. As in the 



calculation of jet cross sections |14j|, the singularity structure of subtraction term is much 



easier to deal with than that of the isolated cross section itself, at the order in perturba- 



tion theory in which we are working. To establish notation, we derive explicit expressions 
for the isolated photon yield in lowest order (O (a° m ) , O («°)). In Section |]TJ], we present 
our derivation of the three one-loop contributions to the subtraction terms. We work in 
n = 4 — 2e dimensions in order to display singularities explicitly. The hard-scattering matrix 
elements are identical to those for the fully inclusive case discussed in Ref. [TLJ, but the inte- 
gration over momentum variables of the unobserved final-state partons is restricted by the 
requirements of photon isolation. All of our calculations are done analytically. In Section IV, 
we summarize our final expressions for the isolated photon cross section E 1 da l e s £ e -^ x /d 3 £, 
and we discuss the logarithmic divergence, £n\(l + e^) — l/sc 7 |, associated with the limit- 
ing case in which x 1 approaches 1/(1 + e^). Numerical estimates, implications for isolated 
prompt photon phenomenology in e + e~ and hadron-hadron reactions, thoughts for further 
work, and suggestions for comparisons with data are collected in Section V. 



II. FACTORIZATION AND LOWEST ORDER CONTRIBUTION 

In this section we drive a general factorized formula for the cross section for isolated pho- 
tons, and we end the section with a computation of the lowest order O {a° em a s ) contribution 
to the isolated energetic photon yield in e + e~ — > 7X. The general notation in this paper is 
the same as that in Ref. [|TJ. 



A. Factorized Form for the Isolated Cross Section 

With the assumption of parton-to-photon factorization, the cross section for an m parton 
final state in e + e~ — > jX, as sketched in Fig. || is 

dPS {m) ■ dzD c ^(z). (1) 

Parton c = 7, <?, <?, and z = E^/E c . The factor dPS^ is the m-parton phase space, and 
D c -+y(z) is a function that describes fragmentation of parton c to the photon. This general 
expression is valid for both inclusive and isolated cross sections. The difference between the 
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inclusive and isolated cross sections resides in the phase space integration for the final state 
partons. For isolated cross sections, because of the isolation condition, not all partons can 
be integrated over all phase space. For example, partons with energy larger than €} l E 1 are 
excluded from the cone of isolation about the observed photon. 

As we proposed in Ref . || , we write the isolated cross section as the following difference 
of cross sections: 

J iso rlrr ind rlrr sub 

7 dn ~ 7 dn 7 dn ■ [ ) 

In Eq. (p|), E^da mcl /d 3 £ is the cross section for inclusive photons. It is well-defined in QCD 
perturbation theory and can be expressed in the following factorized form |l] 



where x c = 1E c j y/s, and the sum extends over c = 7, g, g and g. The short-distance hard- 
scattering cross sections Ecda™^- _> cX / d 3 p c in Eq. (|) are derived in Ref. through one-loop 
level. 

Cross sections for isolated photons E 1 d<T t e s f e -^~ /X /d 3 £ measured, e.g., in LEP experiments, 
are well-behaved and finite. Therefore, the theoretical subtraction term E 1 da^+ b e -^ x /d 3 £, 
defined in Eq. (H), should be well-behaved as well. Since the available phase space for the 
isolated photon cross section is smaller than that for inclusive photons, the cross section for 
isolated photons should be less than the corresponding inclusive cross section, 

E ^ dn - E ^ dn ■ (4) 

Consequently, the subtraction term, E 7 dal+ b e -^ x /d 3 £, should be positive and finite. In 
terms of the definition given in Eq. (Q), E 1 da°+ b e -^ x /d 3 £ can be viewed as a "cross section" 
for a photon "jet" with photon momentum £ and hadronic energy in the "jet" cone E^ ne 
restricted to be larger than E max = e h E^. The virtue of Eq. (0) is that the infrared and 
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collinear singularities of the subtracted term a sub are much easier to deal with, at the order 
in which we are working, than those of a lso itself. 

We assume as a working hypothesis that the subtraction term crl+ b e -^ x can be factored 
in the same way as the inclusive cross section and expressed as a convolution: 

In Eq. (d), we include explicit dependence on the cone size S in the fragmentation functions 
in order to point out that, in principle, the fragmentation functions extracted from the 
cross section for isolated photons may depend on the definition of the isolation cone. Since 
the inclusive cross section a ind in Eq. (0) is well-defined in QCD perturbation theory, it 
is our principal task in this paper to show to what extent the short- distance subtraction 
terms a*+ b e -^ cX are free from infrared and collinear divergences. As we show in Sec. IV, the 
conventional factorization Eq. (|5]) fails in the neighborhood of a well defined point in phase 
space, x 1 = 2E^/ a/s = 1/(1 + e h ). 

The limits of integration over z in Eq. (|5|) are fixed by kinematics: x c = x^/z < 1. 
However, the isolation condition imposes, in addition, a requirement on the total hadronic 
energy in the isolation cone, E^ ne > E mSuX = e^Ey. Generally, E^ ne has two sources: 
Efrag + E partons , as is illustrated in Fig. ||. The fragmentation component, Ef rag , is the 
hadronic component of the "jet" from which the energetic photon itself emerges. In the 
collinear approximation, 

E frag = (1 - z)E c = 0—^) E T (6) 



The second component, E partons , is contributed by other final state partons that are emitted 
into the region of phase space defined by the photon isolation cone. For the subtraction 
term, 

J-'fl Epartons T frag 

Epartons ( ) ^ E mSLX = GfoEy. (7) 



If z < 1/(1 + eh), or, equivalently, Ef rag > ehE^, the constraint Ef l one > ey l E 1 is satisfied 
for any value of E partons . Correspondingly, there is no restriction on the phase space for 



accompanying final state partons. Consequently, a 



sub 



iricl 



for z< 1/(1 + e h ). 



If z > 1/(1 + €h), to satisfy i?™ ne > e^i? 7 , it is necessary that 



Epartons ^ E max E 



frag 



Emin \Z) 



Ey. 



(8) 



We remark that E min (z) > as long as z > 1/(1 + e^), and E min (z) — if z — 1/(1 + e&). 
Taking into account the constraints from the isolation condition, we can rewrite Eq. (| 

as 

da sub 



e+e~— >jX 



E 

c 

+ E 

c 

= E 



z 2 c <i 3 p c 
df 

z 2 c <i 3 p c 



Epartonsl^l ^min 



, 1 )/r:-(--^) 



(9) 



do- 



me/ 



drf sub 



It is important to note in Eq. (||]) that the short-distance subtraction terms E c da^+ b e - ^ cX / d 3 p c 
are needed only for x c < min[x 7 (l + e^), 1] < 1, if x 7 < 1/(1 + e^). 

The modified convolution signs "Cg>" and "<8>" in Eq. (§) are defined in the same way as 
the convolution sign "(g)" in Eq.(|3]), except for the limits of the ^-integration. For functions 
A(x c ) and B(z), we define 

- 1 dz 



A(x c ) ® 5(2) 



A(ar c )®.B(z) = 
A(x c )®B{z) = 



dz 



A(x c = x~/z) B(z) ; 



1 



z 



2 A(x c = Xj/z) B{z) ; 
2 A(a; c = x 7 /z) 5(z) . 



(10a) 
(10b) 

(10c) 



Notice the identity £g> = <g> + <g>. 

Substituting Eqs. @ and flj^) into Eq. (Q), we derive a simplified expression for the 
isolated cross section 
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J ISO 

7 dH 



E 



d 3 p c 

.sub \ 



Ap.sub 
E m e+e-^cX 



^2 I E e+e ~~> cX — E '' ' ~ r v 



d 3 p c 
®D c ^(z) . 



D c ^(z) 



^partons^H Em.in 



d 3 p c ' d 3 p c 

In deriving Eq. (|ll|), we assume for simplicity D^ Jz, 5) = D c ^(z). 



(11) 



When c = 7, D 



c— >7V 



5(1 — z) through order 0(a em ), and E_ 



frag 



(i.e., E n 



€h.E<y). Therefore, we may rewrite Eq. (|Tl|) more explicitly as 



.EL 



e+e — fyX 



E~ 



da 



e+e — fyX 



E s 



dH 7 ^ c= ^_ s d , K 

The isolated short-distance hard-scattering cross sections are defined as 



(12) 



r e+e-^ 7 X 



= ^ K 



AprSub 



d 3 i 



partons inside cone; 
■&partons'>€y l E'y . 



Afriso r1rr incl rlrr suh 
771 a<J e + e -->cX — Tp a(T e+e-^cX rp a<T e+ e --^cX 
t-lr 7T^ = I-Jr ^ tilt 



d 3 p c 



d 3 p c 



d 3 p c 



partons inside cone; 
^partons>E m ^ n . 



(13a) 
(13b) 



The value of E min is specified in Eq. @ . 

Anticipating a result to be derived below, we remark that when x 7 — > 1/(1 + e/J, the 
subtraction term a sub and, consequently, the isolated cross section a lso develop a logarithmic 
singularity, £n\(l + eh) — l/x 7 |. We discuss this singularity in detail in Section [IVD . 



B. Derivation of the Lowest Order Contribution 



In lowest order , O (a° m a°), photon production occurs only through the quark or anti- 
quark fragmentation process sketched in Fig. [|. In this case c = q,q, and x c = 2EJ a/s = 1 
in Eq. ([1]]); and there is no direct production of photons. Owing to momentum balance, the 
quark and anti-quark have equal but opposite momentum in the overall center-of-mass sys- 
tem, and there can be no accompanying parton in the isolation cone around the fragmenting 
parton, i.e., E part 

d& (0)sub 



inside the cone. Therefore, at O (a® m a®), 



E r 



+cX 



d 3 p c 



partons inside cone; 
^partons>E m - in . 



. 



(14) 
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Substituting Eq. (|14"D into Eq. ([□]), we derive 



E.. 



da 



(O)iso 



1 

<l+e h 



dr AtM ind 



c=q,q 



d 3 p c 



X 



(15) 



The lowest order hard-scattering cross section E c da^+ l e n f l ^ cX /d 3 p c in n dimensions is derived 



in our earlier paper [flj : 



, » (O)inci 
E <fo e +e--+cX 

d 3 p c 



-F c PC (s] 
s 



2 / 4^ 



s/4)sm 2 e c ) r(l-e) 



(l + cos 2 # c ) - 2r 



o(x c - 1) 



Xr 



(16) 



with e = (4 — n)/2, and c = q,q. The angle C is the scattering angle of parton c with respect 
to the e + e~ collision axis in the overall center-of-mass frame. The normalization factor is 



[2/s)F Q pc (s) is i 



?PC, 



e 2 + (\v e \ 2 + \a e \ 2 



(1 



|2 1 I |2 



s - M 2 ) 2 + Mf r| 



s( S -MD 



s - M 2 ) 2 + Mfr| 



(17) 



In Eq. (|T7|), both 7 and Z° intermediate state contributions are represented, including 
their interference. The vector (y) and axial-vector (a) couplings are defined through 
ie^iyf + a/75), the vertex coupling between the intermediate vector boson and the ini- 
tial/final fermion pair of flavor / 

At this order, the cross section is manifestly finite in the limit e — ^ 0, and we may set 
e = directly in Eq. flTB"!). Nevertheless, Eq. ( Jl6| ) expressed in n dimensions is valuable for 
later comparison with the higher order cross section. Substituting Eq. flTBD into Eq. (|H>|). 
we obtain the lowest order isolated cross section M 



da 



(O)iso 



d 3 t 



2E 



-F q PC ( 
s H 



a 2 m N c (l + cos 2 0y)— -D 9 ^ 7 (x 7 , /ip), 



(18) 



for x 7 > 1/(1 + 6^). At this order, the isolated cross section vanishes if x 7 < 1/(1 + e h ). The 
angles # 7 and 6 C are identical since we take all products of the fragmentation to be collinear. 
The overall factor of 2 in Eq. ( |T8| ) accounts for the q contribution. Our result in Eq. flTS] ) is 
consistent with those derived in Refs. [EU31. 
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III. SUBTRACTION TERMS AT ONE-LOOP ORDER 



Having calculated the cross section for inclusive photon production in Ref. [|l[], we must 
calculate the subtraction terms defined in Eqs. (Q) and (|1^) in order to derive the com- 
plete expression for the isolated photon cross section. As in the inclusive case, there are 
three distinct contributions to the subtraction terms for e + e~ —>■ 7X at one-loop level in 
perturbation theory: 

e + e - _> 7; 0(a em ) (19a) 

e + e~ — > q (or q) — > 7, 0(a s ) (19b) 

and 

e+e" g -+ 7. 0(a s ) (19c) 

In Eqs. ( |19b|) and ( |19c|) , we have in mind contributions from quark (or antiquark) and 
gluon fragmentation to photons in the three-parton final state process e + e _ — > qqg. The 
first contribution, Eq. ( |19a| ), arises from e + e~ — ► qq-f where the 7 is not collinear with either 
q or q. 

In this section we derive and present explicit expressions for the contributions to 
Er r dal+ b e -^ lX /d 3 £ from each of the three processes in Eq. ([T9|). 

A. Factorized Form for the Short-Distance Hard Parts: a sub 

To calculate the short-distance partonic cross sections (hard parts), a sub in Eq. (|13|). 
we first apply the factorized form in Eq. (^]) to parton states perturbatively order-by-order 
in the coupling constants, and we then extract the perturbative expressions for a sub . We 
illustrate this procedure in this subsection for each of the processes in Eq. ( f[9|) 

The Feynman graphs for c^+J^ x are sketched in Fig. [|. To derive the direct production 
term corresponding to Eq. ( |19a| ), o^+e-^yX' we apply Eq. (|9]) perturbatively to first order 
in a em . We obtain 
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(1)«*6 _ a . 



^ 7 x(^)^%(z) + ^^(x^D^z) 



E q (orE<i)>e h E-, 



+ (g - g) . (20) 

The zeroth order subtraction term 5"e+e--»cX vanishes, as shown in Eq. JT^). Since 
the zeroth order photon-photon fragmentation function Dj® 1 {z) = 6(1 — z), and 
^.(i)mc^ ^^x)^ Jz) vanishes, the expression for the short-distance hard part becomes 



a e+e-^jX\ X j) — "e+e-^X 



where the terms with the minus sign are often called collinear counter-terms. In this equa- 
tion, o'f+^L.gx is obtained from Eq. ([U]), and the modified convolution "®" is defined in 
Eq. ( |10c| ). The perturbative fragmentation function D q l \ (z) in n-dimensions can be calcu- 
lated by evaluating the diagram in Fig. |6|. We obtain 

and d£>(z) = Dg)(z). Although 

ff e + e -_, 7 j( and the perturbative fragmentation functions 



Dg? OT g\_^ 7 (z) are both formally divergent as e — ► 0, these divergences cancel and leave a finite 



expression for d'^^ b x , if the conventional QCD factorization theorem holds ||T0| . 



For the gluon fragmentation contribution: e + e _ — > g — > 7 (Eq. (|19cf) ), we apply Eq. (g) 
to a gluon state perturbatively to first order in a s . Letting E g be the gluon's energy, we 
obtain 



(l)«d> 
a e+e-->gX 



+ ^ q xi^)®Dfl g {z) + a^xix^D^z) 

+ (<Z - ?) , (23) 

where x g = 2E g /y/s, and x q = 2E q /^fs. In Eq. (^3|), the modified convolutions are defined 
in the same way as in Eq. fllOD, except that x 7 is replaced by x g , and is replaced by e m i n ; 

/ % -^min(^) 

£min{Z) = — 

= (l + e h )2-l>0, (24) 
13 



where z = x^/xg. Because z < 1, we note the restrictions 

emin(^) < t h ; or — > — — . (25) 

-L + e m in J- + £h 

Since D g °X g {z) = 5(1 — z), and &e+^ b -+ q x( x q) = 0, Eq. (|23|) can be rewritten as 



e+e-^gXV X 9J ~ U e+e~^gX 



(26) 



In Eq. (^6|), the divergent cross section cr^^- b x is evaluated from the Feynman graphs 
shown in Fig. [7|, and the quark-to-gluon collinear divergences are embedded in the first- 
order fragmentation function D^ g . The function D^2, g is the same as D£) of Eq. (p2| ), 
except that e;j(a: em /27r) is replaced by Cj?(a s /27r). The color factor Cp = 4/3. 

For the quark fragmentation contribution: e + e~ — ■> g — ► 7 (Eq. (|19b|) ), we apply Eq. (^|) 
to a quark state perturbatively at first order in a s . In a fashion similar to the derivation of 
^e+e--*gXi we derive the short- distance hard part for quark fragmentation 

The Feynman graphs that contribute to cr^^ b x are sketched in Fig. || In this fragmenta- 
tion process, the quark is effectively "observed" through q — > 7 fragmentation. For the real 
gluon emission diagrams, sketched in Fig. [8]a, the gluon and antiquark are not observed and 
their momenta will be integrated over. Since e min (2;), defined in Eq. (^4[), can be equal to 
zero, the contribution of the real emission diagrams to 0^+^ x shows both infrared and 
collinear singularities. The infrared singularity associated with soft gluon emission should 
be canceled by a contribution from the virtual gluon exchange diagrams shown in Fig. §b, 
if the conventional QCD factorization theorem holds. When the gluon is parallel to the 
fragmenting quark, the real emission diagrams manifest a collinear singularity that should 
be canceled by the negative term in Eq. fl2"7|). 

The perturbative fragmentation function DjQJz) in Eq. ( p7|) is determined from Feyn- 
man diagrams sketched in Fig. |9]. In n- dimensions, we obtain ]TJ 

1 



14 



-e 



(28) 



where the "+" description is defined as usual, 



rb) + s rb- , < 1 -'>r <fa 'rb- (29) 



B. Parton Level Cross Sections: a^ sub 

In order to derive the short- distance hard parts at one- loop level, defined in Eqs. (pT|). 
(p6|) and fl2"T|), we must compute the formally divergent partonic cross sections "g+^-^cX m 
n-dimensions, for c — j,g,q and q. We must evaluate Feynman diagrams for e + e~ to three 
particle final states, (qqj) or (qqg)] and e + e~ to two particle final states (qq) with one-loop 
virtual gluon exchange. As derived in Ref. Q, the general expression for the partonic cross 
section is 



da (i)sub = y \If? c (s)] e 2 C q \(H X + H 2 ) dPS {m)su \ (30) 
— ' s 4 

where m = 2 or 3 corresponding to the number of final state particles, and the constant C q 
is an overall color factor. In Eq. fl30|), the functions H\ and i?2 are defined by 

H x = - 9lw H", H 2 = -hp. H^. (31) 

Function H^ v is the hadronic tensor, and q^ (/c M ) is the sum (difference) of the incoming 
e + and e~ momenta. The factor dps^> sub in Eq. ( pOf ) is the multi-particle phase space 
element. In n = 4 — 2e dimensions we obtain JTJ 

,po(2M 1 1 rf W 4 ^ V 1 27r ^(a;, - 1) 

2(2tt) 3 ^ ^(s/4)sin 2 0j r(l - e) s z< 1 J 

for the two-body final state with parton i being observed; and 



dPS {3) 



sub = 1 1 d s Pi ( 4vr V 1 2n 5 ( Xi - (1 - y jh )) 
"2(2tt) 3 V(s/4)sin 2 ^J r(l - e) s x { 



s 

X 4 



lYf^Y 1 ^n- 3 fe) 



2tt/ V s / rfl - e) a,- 



n— 3 



x — 1 — 0(1- ?/i2 - yi3 - 2/23) , (33) 

2/12 I/13 I/23 
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for the three-body final states. In Eq. (j33|), we use labels "1", "2", and "3" for the final- 
state q, q and 7 (or g), respectively; and we use label "i" (=1,2 or 3) to designate the 
"observed" one. Symbol "j" (7^ i) labels the particle over whose momentum we integrate, 
and the momentum of i 7^ j) is fixed by the momentum conservation 5-function. In 

Eq. (|3"3"|), fi„_3 = 27r (1 ~ <E )/T(l — e). The dimensionless invariants %i and in Eq. ( |3"3"D are 
defined as 



2Ei 



X j 



Vij 



2Pi ■ Pj 



(34) 



with i, j = 1, 2, 3. 

The multi-particle phase space expressions in Eqs. ([32] ) and (|33"D, are formally identical to 
those for our calculation of the inclusive cross section QXJ , except that the limits of integration 
differ here, owing to the isolation condition. These phase space limits are derived below for 
each individual subprocess. 

As indicated in Figs. ^, [7] and |8|, all three one-loop partonic cross sections have the same 
e + e~ to three-body matrix elements. However, the quark fragmentation cross section has an 
extra contribution from virtual gluon exchange diagrams. We derive the functions Hi and 
H 2 appearing in Eq. (^) by evaluating the square of the matrix element for the Feynman 
diagrams in Fig. [5] (or Fig. 0, or Fig. ||a). 



i2l = .S(l-r){(l-r 
H 2 



-4|(l-e) 
4 



2/13 | Z/23 

2/23 2/13 
2/13 | 2/23 
2/23 2/13 



2/132/23 



2,2 
Vlk + 2/2/c 



+ 2 
4e 



+ 2 
2/12 

.2/132/23 
Vlk- 



2/12 
2/132/23 



— e 



(35) 



2/132/23 

We temporarily omit overall coupling factors of (e//) 4 for (qq'j) and (e//) 2 (g/i e ) 2 for (ggg). 
The invariants j/jfc with 2 = 1, 2, 3 are defined as 

2pi ■ k 



(36) 



The identity 
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yik + V2k + 2/3fc = o 



(37) 



follows from the equalities k ■ q = and q = p\ + p 2 + p$. 

We choose to work in the overall center of mass frame, sketched in Fig. [0], when integrat- 
ing Hi and H 2 over three-body phase space. Letting 9 ik be the polar angle of the observed 
parton momentum pi with respect to the momentum k, and angle 9j X the n-dimensional 
generalization of the three-dimensional azimuthal angle 0, defined through pj, we derive 
general expressions for the three pfl. 



Vik = -Xi cos9 ik ; 
Vjk = 

Vhk = 



'VihVjh - 




cos 6 ik - 




Xi 








'VijVjh - 


Vih 


cos 6 ik + 




Xi 






Xi 



sin 6 ik cos 9a 



sin 9i k cos 9j X . 



(38) 



These general expressions satisfy the identity in Eq. (37). In the reference frame we are 
using, 



dttnsipj) = dcos9 jx (1 - cos 2 9 jx ) n 2 rffi n _ 4 (p j ) . 



(39) 



Since we choose the photon direction to define the z-axis ( Fig. |10|), with our cone 
definition, the range of the azimuthal angle 9j x integration should not be affected by the 
isolation cut. Therefore, just as the inclusive case [|IJ, the integral over dcos9j x is done from 
cos9j X = — 1 to +1 when H 2 is integrated over phase space. The expression for three-body 
phase space, Eqs. (^) and (|39|) , is an even function of cos^. Correspondingly, terms in 
H 2 that are odd functions of cos 9j X do not survive after integration. Because H 2 (Eq. (|35|) ) 
depends only on the square of y\ k and y 2k , after eliminating all terms linear in cos 9j X , we find 
that the only remaining dependence on 9j X in H 2 is of the form cos 2 9j X . We may integrate 
over 9j x independently of other variables, or, equivalently, we may replace cos 2 9j X in H 2 by 
its average in n-dimensions and thereby eliminate 9j x dependence in H 2 completely. 



The average of cos 2 9j X in n-dimensions is pi] 



<C COS 9j x i> 7i— dim" 



2(1 



(40) 
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We use Eq. (f40|) to replace cos 2 9j X in H 2 by 1/2(1— e), and we denote the resulting expression 
HI". 

In following subsections, for a given observed particle (e.g. 7 or q) , we derive expressions 
for yf k , for i = 1,2,3, substitute these yf k into Eq. (|35|), replace cos 2 ^ by 1/2(1 — e) to 
obtain (Hi + H^"), and finally, substitute these (H1 + H2 ) into Eq. fl3U| ) to get the partonic 
cross sections: a^ sub . 



C. Derivation of a^ l ) svh v 



In this subsection we present an explicit derivation of the finite hard-scattering cross 



section a 



(l)sub 
e+e — — >yX 



at order a em , as expressed in Eq. (|2T|). 



Since the photon is the observed particle, we set i = 3. The assignment of j and h is 
arbitrary and has no effect on the final result owing to symmetry between the quark and 
antiquark. We choose j = 1 and h = 2; i.e., "1" for the quark, and "2" for the antiquark. 
Referring to Eq. (|34l) , we remark that j/13 is proportional to the invariant mass of the photon 
and quark system. Starting from Eq. (|38D , and after using Eq. (|40D, we obtain the effective 
expressions 



2 2 2/1 



Vik 



1)\k 



2/232/12 - 2/13 



^3 

2/132/12 - 2/23 



COS 2 #3 + 



x 3 



cos 2 6» 3 + 



1 - e 
1 

1 - e 



2(2/122/132/23; 
3 

2(2/122/132/23) 
£3 



Substituting these y 2 fc 5 with « = 1, 2, 3, into Eq. (p5|) , we derive 



7 (#i + # 2 e// ) = ^(e//) 4 ( (l + cos 2 9 3 - 2e 



sin 2 0, 



sin 2 0,. 



V2/23 2/13 



M 



(41) 



+ (1 -3 cos 2 ^ 3 



42/ 
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eg(e/i e ) 4 j (l + cos 2 # 3 -2e 
+ (l + cos 2 3 - 2e 



l + (l-x 3 ) 5 



(-■ 




V2/13 


2/23/ . 



2/12 
2/132/23 



K-- 




V2/13 


2/23/ 



(l - 3 cos 2 fl 3 



4(1 -x 3 ) 



4 



(42) 



In Eq. (H), we introduce the overall coupling factor and neglect terms that do not contribute 
in the limit e — > 0. The variable yj 3 with j — 1, 2 is defined as 7/7-3 = yj 3 /x 3 ; x 3 is the same 
as x 7 , and 9 3 is equal to # 7 in this case. 

The two 6 functions in dPS^ 3 \ Eq. fl33|) , are used to do the dyn and g?7/ 2 3 integrations, 
and 7/13 is left as the integration variable. Different choices of the integration variable are of 
course equivalent. In terms of 7/13 and x 3 , we have identities 

7/12 = 1 - £ 3 , 7/23 = x 3 - 7/13. (43) 

The limits of integration over y 13 are derived from the isolation requirement for the 
subtraction terms. In our parton level approach, hadronic energy in the isolation cone 
means a parton must be present inside the isolation cone of observed photon. Second, this 
hadronic energy must be larger than 6hE y . We use 5 to denote the half angle of the isolation 
cone. The statement that the quark (parton "1") is inside the isolation cone provides the 
relationship 

„ cos 2 (V2) 

2/23 > ~ ■ 2 (X/oV 44 

1 — x 3 sin (0/2) 

whereas the statement that the antiquark (parton "2") is inside the isolation cone leads to 

. cos 2 (V2) 

2/13 > ■ 2/r/ V 45 

1 — x 3 sin (6/2) 

Using Eq. ( |4*3] ) and < 7/13 < x 3 , we derive the condition that either the quark or the 
antiquark is inside the isolation cone: 

(l-x 3 )sin 2 (V2) 

< 2/13 < "j ■ 2fX ,o\ ; 

1 — x 3 sm (6/2) 

cos 2 (5/2) A 

< 2/13 < 1. (46) 



1 -a; 3 sin 2 (5/2) 

These two regions of the 7/13 integration do not overlap as long as x 3 > 1 — cot 2 (6/2). 

The second part of the isolation constraint is the requirement that the parton's energy 
inside the isolation cone be larger than e^E 1 . We derive 
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r,w, w m . (1 -x 3 )sin 2 (o72) 

maxO, 1 + e,- l/x 3 < y 13 < ^ . 2 ^/ 9 / 

1 — x 3 sm (o/2) 

, (l-x 3 )sin 2 (5/2) 



< y 13 < min[l, 1 - (1 + e h - l/x 3 )} . (47) 



1 -a; 3 sin 2 (5/2) 

For simplicity of notation and to facilitate comparison of our results with those of Ref . || , 
we define 

1 — x 3 sm (0/2) 4 
y m = l + e h . (48) 

Henceforth in this subsection, we set x 3 = x T The first of the conditions in Eq. (|47|) indicates 
that y m should be less than y c . Therefore, the subtraction term should vanish if y m > y c , 
i.e., 

*(1)«*6 

£ 7 e+ ^ 3 ^ 7X =0 if x 7 > x™ x (5, e h ) . (49) 

Here, x™ ax (<5, e^) is expressed as 

x- ax (5, Cfc ) = (— |-) 2 9 (50) 

V 1 + e ft 7 ! + A _ Aeh siri 2 (5/2 )/(i + e/i )2 



We note that 



if 5 ^ 0. We define 



< ax (<W)>-^ (51) 



A = x™ (5, Cfc ) - — L- ~ e h j , (52) 



where we expand x™ ax (<5, e/J of Eq. ( pT] ) for small 5. For 5 = 20°, typical of LEP experiments, 
A ~ O.OSe^. The interval between the critical value x 7 = 1/(1 + 6^) and x™ ax is consequently 
very small. Equation ( f4D|) is actually a condition due to energy-momentum conservation. 

We consider next, in turn, the regions x 7 < 1/(1 + e/J and x 7 > 1/(1 + e/J. If x 7 < 
1/(1 + e^), the integration over y 13 for (icr^^ su6 has two separate intervals, defined in Eq. 
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/13 



He 



dyi3 + / 
Ji- 



/13- 



(53) 



10 Ji-yc 

Substituting Eq. (|42"D into Eq. (^3), and performing this j/13 integration, we derive parton 
level cross section 



da 



(l)sub 

~dH 



2E 



F q PC (s) 



2 AT 



47T/i 2 



1 



x <^ (l + cos^ 7 -2e) 



l + (l-s 7 ) 2> 



+ (1 + cos 2 7 ) 



+ (1 - 3 cos 2 ^ 



1 + (1 - x^f 

+ x 7 (l - 2y c 



In 



^/4)sin 2 # 7 / T(l-e 

1 

e 
s 



1 



2 / ^em 
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+ £n(x (1 - x 7 )) + £n(y c ) + y c 



MS' 



(54) 



In Eq. (|54|), the usual modified minimal subtraction scale is 



-IE 



(55) 



where 7^ is Euler's constant. The 1/e poles in Eq. (|54)) arise from the 1/^13 and 1/^23 terms 
in Eq. P2]). Referring to Eq. (|2"T|) , and using Eqs. (|TB|) and (|2"2]), we observe that the 1/e pole 
in Eq. (Q) is exactly canceled. Consequently, for x 7 < 1/(1 + e^), the finite short-distance 
hard part is 



da 



(l)sub 
e+e~— >-yX 

~dH 



2E 



a 1 N - 



X 



7 J 



a. 



2tt 



x (1 + cos 2 6> 7 



1 + (1 - 



in 



+£n (1 - x 7 ) 



5 1 



+ in(xl{\ - x^)) 



(56) 



In Eq. (|56|), we neglect terms of 0(5 2 ). 

If x 7 > 1/(1 + e/i), the negative term, defined through the (g> convolution in Eq. (PT|). 
vanishes. In this region of phase space, there should not be any collinear divergences in 
the partonic cross section, E 1 da^+^,- b ^ x /d 3 £, since there is no counter term to cancel them. 
If the fragmentation process were exactly collinear, the subtraction term E 1 da^™ b x / d s £ 
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would vanish, kinematically, for x 7 > 1/(1 + e^). However, a finite isolation cone 8^0 
allows for a non-vanishing E^dcr^^* x /d?£ even when x 7 > 1/(1 + e/J. Equation (f49|) 
shows that there is a narrow interval, 1/(1 + 6^) < x 1 < x™ ax (5, e^), in which a non- vanishing 
Eryda^l-^^x I d 3 £ is allowed kinematically. 

In the narrow interval, 1/(1 + e/J < x 7 < x™ ax (<5, e/J, the integration over y 13 for da^ sub 
has two separate regions: 



dy 



13 



dy 



13 



Vm 



1-3/c 



#13 • 



(57) 



Equation (f48|) shows that y m > 0. Therefore, $13 > and f/13 < 1. Consequently, there 
is no collinear divergence and the integration over y 13 can be done completely in n = 4 
dimensions. 

For 1/(1 + th) < Xj < x™ ax (<5, eh), we derive 

, (l)sub 



E, 



' dH 



2E 



x 1(1 + cos 2 6> 7 ) 



+ (1 -3 cos 2 ft 



a 2 N — 

em x 



1 + (1 -x 7 ) 

A x i 

- 4x 7 (y c - y n 



2tt 



2E 



a 2 N 



£n— + £n 

Vm 



(Vc - Vm) 



1 - Vrr. 
1 - Vc 



2tt 



x(l + cos 2 6> 7 
+0(5 2 ). 



1 + (1 - x n 



6z 



(1 - x 7 )£ 2 /4 
l + e h - l/x~ 



(58) 



We remark that E 1 da^- b _^ yX /d 3 £ in Eq. (|58D manifests a logarithmic divergence, of 
the nature of a collinear divergence, as x 7 — > 1/(1 + e^). This problem arises from the 
incompatibility of collinear fragmentation, used to define the counter term in Eq. fl21~|), and 
the cone fragmentation used in the definition of the partonic cross section E 1 da^™ b x / d s £ 



in Eq. fl2T|). More discussion of this issue is presented below, in Sec. [IV C 
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D. Derivation of a 



(l)sub 
e+e~ —+gX 



In this section we present our explicit expression for the finite short-distance hard part 
for the order a s gluon fragmentation process, e + e~ — > g — > 7, expressed in Eq . (^6|) . The 
similarity between the Feynman diagrams for e + e~ — > 7, shown in Fig. |5], and those for 
e + e~ — > g — > 7, shown in Fig. [7] allows us to exploit the results derived in the previous 
subsection. 

We derive a[+H- b ,„ v from the expression for &Q™ b x , given in Eqs. fl56|) and fl5"8|), by 



making the following four replacements: x 7 — > x g ; N c — > N c Cp] e 2 e 2 of the final photon 
emission vertex by g 2 = Arrets] and — > e m j Jl (z). The last replacement is absent for the 
inclusive case. Here z = x 1 /x g . 

If x g < 1/(1 + e min ), which is the same as a; 7 < 1/(1 + these replacements in Eq. 
provide the finite short- distance hard part 

. (l)sub 



E, 



d K+e-^gX 

d 3 Pg 



2E 



2 AT 



X 



91 



X (1 + cos 2 a 



in 



+ in{x 2 {l-x g )) 



+tn (1 — x 



+ x 



+0(S 2 ) . 



(59) 



If x g > 1/ (1 + e min ), which is the same as x 7 > 1/(1 + e/J, we make the four replacements 
in Eq. (|58|) and find 



E n 



da 



(l)sub 
e+e~-*gX 



d 3 Pg 



2E 



q 



x(l + cos 2 63) 
+0(S 2 ). 



a 2 N — 

X 9 1 



1 + (1-Xgf 



(l-x^ 2 /4 



(1 + e h )(xj/x g ) - 1/x 



9/ 1 



(60) 



From energy-momentum conservation, we derive an equation similar to Eq. (E9| 



da 



(l)sub 
e+e--^gX 



d 3 Pg 



if x 9 >x max (z,5,e,). 



(61) 
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The maximum value x max is 



x max (z, 5, e h ) = — ; . (62) 

A 1 + e h) 1 i ,A| 4 S in 2 (V2)( 2 (l+e h )-l) 

The gluon fragmentation contribution to the 0(a s ) subtraction term in o^+g-^, ~ is 



provided by the convolution 



(l)snfc 

1 ' z 



7 (l)suo 



max 



(l)sub 

:+e~- 

r/ :i /' 7 V 9 z 



d 3 £ 

If collinear fragmentation is assumed for g — > 7, #3 = # 7 . We remark that, after convolution 
with the gluon-to-photon fragmentation function, the subtraction term E g da^+^- b _^ gX /d 3 p g 



in Eq. (pO|) develops a logarithmic divergence £n(l + e/i — as a; 7 — ► 1/(1 + £/,). More 

discussion of this issue is presented below. 



E. Derivation of q-'-V™ 6 



qX 



In this section we evaluate the finite short- distance hard part ^e+e-^qX a ^ or der a s , 
defined in Eq. (f2~T|), for the quark fragmentation process e + e~ — ► g — ► 7. Because z > 
max[x 7 ,l/(l + eft)] in the convolutions of Eq. (|Tl|) , our general analysis of Sec. |T| shows 
that the short- distance hard part o , ^+^ h x is needed only for x\ = x c < min[x 7 (l + eh), 1]. 
We recall that x\ = 2E q / ' y/s. We analyze, in turn, the intervals x 7 < 1/(1 + e h ) and 



Xj > 1/(1 + Eh). Discussion of the special case x 1 = 1/(1 + e/J is reserved for Section [IV D 
If Xj < 1/(1 + €h), then x\ < 1 for all values of z in the convolution. In this interval, the 
virtual diagrams of Fig. ||b, whose contribution is proportional to 8(1— x±), do not contribute 
to ffg+g-Lx- The sc l uare of t ne matrix element for o , ^^ b x is therefore the same as that 
in Eq. fl35|), except that the yf k variables are no longer those of Eq. (f4"T|). Because the quark 
is now the fragmenting parton, the labels are instead i — 1, j — 3, and h = 2. We obtain 
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2 2 2 /1 

Vik = x i cos 0i ; 



2/2fc 



2/3fc 



2/132/23 ~ Vl2 
X\ 

2/122/23 ~ //[:; 



cos 2 6*i 



1 



1 2 



cos 2 #i 



1-e 
1 



2(2/122/132/23] 



-I 



2(2/122/132/23) 

^y»2 
J, 1 



sin 2 6>i ; 
sin 2 6 1 . 



(64) 



In Eq. (|64"D, $i is the angle between the quark momentum pi and the momentum /c. In 
deriving Eq. ([64]), we dropped terms linear in cos8j X , and replaced cos 2 0^ by its average 
in n-dimensions, 1/2(1 — e). Substituting these yf k into Eq. (|35"D, and using the identities 
y 23 = 1 - xi and j/ 12 = x x - 2/13, we derive 

2/12 



e// 
2 



(e//) 2 (<7/^ 2 



+ 



(e//) 2 (<7/^ 



+ 



1 + cos^ 61 - 2c 



2/13 | 2/23 
v 2/23 2/13 



2/132/23 



1 - 3 cos 2 e 1 
1 + cos 2 0! - 2e 
1 + cos 2 0! - 2e 



22/ 



12 



2/13 



i+i 1 + 

1 - Xi) 2/13 1 - Xi 

2 fl — xi 



2/13 



1 — Xi 



2/13 



Xi 



1 - 3 cos 2 0i 



xi V xi 



(65) 



The overall coupling constant (e/i e ) 2 (g//) 2 is restored in Eq. ( |65|) . Since xi < 1 for all values 
of z, there is no infrared divergence associated with 1/(1 — xi) terms in Eq. (|65|). 

The isolation condition for the partonic cross section cr+e-^x requires that either the 
gluon or the antiquark be in the isolation cone of the quark and have energy larger than 
E m[n defined in Eq. (0). Following an analysis similar to the one that led to Eq. (W7h, we 
derive the limits of the integration over 2/13 = y\zjx,\. 



m&x[0,y m ] < yi 3 < y c ; 

1 - y~c < 2/13 < min[l, 1 - y r , 



(66) 



In Eq. (|66|) , y c and y m are defined as 

(1 - x x ) sin 2 (<V2) 
1 - x lS in 2 (5/2) 



Ve 

ym 



(l + e h )z ; 

Xi 



(67) 
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with z = x 1 jx\. In analogy to Eq. fl6"T|), 



da 



(l)sub 



d 3 pi 



if x x > x max (z,S,e h ) 



(68) 



where x max (z, 5, eh) is defined in Eq. (|62|). 

In the region x 7 < 1/(1 + 6^), y m < 0. The integration over ?/ 13 has two separate intervals: 



r rye rl 

/ dy 13 => / dy 13 + / dy 13 . 

J JO Jl-Vc 



(69) 



Integrating over yi 3 , we obtain 



Ey 



da 



(l)sub 
e+e — — *qX 



d 3 pi 



s H 



o? N ( ^ 



1 



x <^ (1 + cos 2 0i -2e) 



^(s/^sin 2 ^; r(l-e) 

\ — X\)\ c 



-c F (^ 

x\ \2ir 



+ (l + cos 2 Ox) 



1 + x\ 

1 — Xi 

+ a-xx) 



in 



J'- 



( b 2 

+ £n(x 2 (l - xx)) +£n[(l- xx) — 



MS- 



(70) 



Terms of 0(5 2 ) are neglected. The 1/e pole in Eq. ( |70"1) arises from the l/yx3 term in Eq. (|65|). 
corresponding to the collinear singularity when the gluon is parallel to the fragmenting quark. 
Combining Eqs. ( [TBI) an d Eq. (|28|) , and using the fact that xx < 1, one may verify that this 
collinear singularity is canceled exactly by the subtraction term in Eq. (|27|). 

For the finite short- distance subtraction term in the region x 1 < 1/(1 + e^), we obtain 



Ex- 



da 



(X)sub 
e+e~-*qX 



d 3 pi 



-F q PC (s) 
s y 



x (1 + cos 2 Bx 



a 2 N — 

"em 1 v c 

X\ 

1-11. 



a 
2n 



CM — 



hi 



+ £n(x\(l - xx)) 



+£n (1 - xx 



S 2 



+ a-xx) 



(71) 



Turning to the case X\ > 1/(1 + e min ( 2 )), which is the same as x 7 > 1/(1 + e/J, we note 
that y m > 0. The integration region over 7/13 = yxz/xi now has the form 



/13 



dy 



13 



Urn 



/13- 



(72) 
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Equation (|68|) shows that x\ < 1 for a nonvanishing Eida^^ b _^ qX / d 3 pi. Since xi < 1 
and yi3 > f/ m > 0, there is neither an infrared nor a collinear divergence in this region. We 
perform the integration of Eq. flB5|) over y 13 in n = 4 dimensions, and we obtain 



da 



(l)sub 
e+e--*qX 



d 3 pi 



a 1 N 



2^ 



x(l + cos 2 9i) 



l + xl 

1 — Xi 



hi 



(l-x 1 )5 2 /4 



(1 + e^^Xj/xt) - l/xi i 



(73) 



Terms of 0(5 2 ) are dropped. In deriving Eq. (|73|) , we used the fact that the second term in 
Eq. ( f27"f ) vanishes in this region. 

To summarize the derivation of this subsection, we present a general form for the 0(a s ) 



subtraction term to a 



da 



(l)swfe 



(l)sub 
e+e~ — >7 X 

1 



dH 



E 



via quark fragmentation 

da { }J™ b 



dz 

z 



E l - 



qX 



d 3 pi 



X\ = — 

z 



D 



• (74) 



For collinear fragmentation q — > 7, ^3 = 9^. 

We remark again, as we did in our discussions of the direct and gluon fragmentation 
contributions, that after convolution with the quark-to-photon fragmentation function, the 
subtraction term in Eq. ( |T4] ) manifests a logarithmic divergence of the form £n(l + — 1/ x 7 ) 
when x 1 — > 1/(1 + e^). 



IV. ONE-LOOP CONTRIBUTIONS TO THE ISOLATED CROSS SECTION 

Combining the one-loop subtraction terms calculated in last section and the one-loop 
contributions to the inclusive cross section derived in Ref. [[!]], we present in this section the 
complete one-loop contributions to the cross section for isolated photons in e + e~ collisions. 

For the convenience of later discussion, we first list all one-loop contributions to the 
inclusive cross section, and then we present analytical results for the isolated cross section 
in three subsections corresponding to x 7 < 1/(1 + e^), x 7 > 1/(1 + e^), and x 1 = 1/(1 + 
€h), respectively. We conclude this section with a discussion of the logarithmic divergence 
associated with the special point a; 7 ^ 1/(1 + e h ). 
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A. One-Loop Contributions to the Inclusive Cross Section 



Complete one-loop contributions to the inclusive cross section, E^da^-^x/cPi, are 
derived in Ref. []]]. The results are summarized as follows. 
One-loop direct production of e + e~ — > 7 at 0(a em ): 



-E-v 



da 



(l)incl 
e+e~—*-yX 

~dH 



2E 



-F q PC (s) 
s H 



x m + cos 2 6> 



* AT 



1 + (1 - X, 



2 ^ Q^em 

2tT 



in 



+ ^n(x 2 (l-x 7 )) 



+ (1 - 3 cos 2 6- 



2(1 -x. 



(75) 



One-loop fragmentation contribution via a gluon: e + e — > g> — > 7 at 0(a s ): 



.EL 



(l)inci 

dH 



1 



(l)iraci 
e+e — — >gX 



(76) 



where the short-distance hard part is 



<i<7 



(l)iraci 



2E 



2 IT 



3 1, 



X (1 + COS 



l + (l-a; 9 ) 2 ~ 



+ Mx 2 (l-x 5 )) 



MS- 



+ (1 -3 cos 2 7 ) 



2(1 



(77) 



with x 9 = 2E g /y/s, and 6> g = 7 . 

One-loop fragmentation contribution through a quark, e + e" —>■ q —>■ 7 at 0(a s ). 



E. 



da% 



(l)incl 

e+ e-^qX^-yX 



E 



1 dz 



j ~ (l)incZ 



d 3 pi 



(78) 



where the short- distance hard part is 



d& (l)ind 



e+e-^qX 



d 3 pi 



'-F„ PC ( 



x 1(1 + cos 2 ft, 



a 2 N - 



2tt 



J^L + 3, (1 _ Xl) ) fe (_f 



MS- 



2 U-zi 
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+6(1 -xx) 



2ir 2 9\ 1 



(3xi - 5) 



+ (l - 3 cos 2 7 ) | . 



(79) 



Here 6\ = 9 1 is used, and the "+" description is defined in Eq. (|29] 



B. One-Loop Contribution to the Isolated Cross Section when x 7 < 1/(1 + e/J 



Substituting Eqs. ( ^6|) and (|75|) into Eq. ( |13a|) for c = 7, we derive the 0(a em ) one-loop 
direct production of isolated photons via e + e~ — > 7: 



e+ e~ —>-yX 



2E 



-Fr(s, 



a 2 N - 



- J 



2 / "em 

2tT 



x |(1 + cos 2 6> 7 
+ (1 -3 cos 2 7 ) 



1 + (1 -x^f 



in 



1 



(1 -x 7 )5 2 /4 / 



2(1 -x. 



(80) 



where we have dropped terms of 0(S 2 ). By integrating over 7 , we may verify that Eq. (|80|) 
is consistent with the analogous expression derived in Ref. @. 

From Eqs. (|TTD, (|59|) and (|77|), we obtain the 0(a s ) one- loop gluonic fragmentation 
contribution to the cross section for isolated photons: 



.EL 



(l)iso 

dH 



dz 



where the short- distance hard part is 



j a (l)iso 

^ ^ r 9 " Z 



D g ^(z,n 2 F ) 



(81) 



(l)iso 



2E 



X 



.'/J 



2^ 



x <Ml + cos 2 



+ (1 -3 cos 2 7 ) 



2(1 -x„ 



(l-2 9 )<5 2 /4 



- x r , 



(82) 



Similarly, we derive the 0(a s ) one-loop quark fragmentation contribution to the cross 
section for isolated photons: 
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da 



(l)iso 

~dft 



E 



dz 



where the short-distance hard part is 



E, 



(l)iso 
e+e-^qX 



d 3 Pq 



-F„ PC ( 



a 1 N - 



j a (l)iso 
E d K+e-^X 



2^ 



D 



>,/4) 



(83) 



| (i + cos 2 e. 

+ fl - 3 cos 2 0. 



1 + x\ 

1 — X\ 



hi 



(1 -x^/A 



3 

2 VI 



{xi - 3) 



(84) 



Equation (|4]) is derived by substituting Eqs. (|7l~l) and ([n|) into Eq. ( |13b|) . We observe after 
integrating over z, that Eq. (|83"D develops a logarithmic divergence £n(l/x 7 — (1 + e^)) as 



1/(1 + eft). This divergence is caused by the 1/(1 — x\) terms in Eq. (|S4]). A detailed 



discussion of this divergence is presented later in subsection [IV D 



Assuming D, 



g->7V 



Dg^(z), the one-loop antiquark fragmentation contribution to the 



cross section for isolated photons is the same as that of the quark contribution in Eq. (|33] 
except that the is over Q- 



C. One-Loop Contribution to the Isolated Cross Section when x 7 > 1/(1 + e/j) 

As pointed out in Section [[II], the subtraction terms a^ sub would vanish for x 7 > 1/(1 + 
e/j) if all fragmentation processes were exactly collinear. But, due to the finite cone size, 
there is a small region of phase space where the subtraction terms are finite. Since the value 
of €h is very small in most experiments, very limited data from e + e~ — > 7X are available in 
this region. Nevertheless, the isolated cross section in this region is interesting theoretically 
and important for understanding the isolated cross section in hadron-hadron collisions. 

Kinematics require the subtraction term E^da^+^ b x /d 3 £ to vanish if x 7 > x 7 nax (<5, e^), 



defined in Eq. ([50"D; therefore 

, ~ (l)iso 

da,l 



E~ 



■yX 



E~, 



da 



(l)incl 
e+e~-»7X 



if x 7 >x 7 nax (5,e,). 



dH 1 dH 

However, if 1/(1 + e h ) < x 7 < x^ 3X (S, e^), we use Eqs. ( |58| ) and (|7|) to derive 
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(85) 



K 



dH 



2E 



1 + COS 0, 



a 2 JV — 
1 + (1 - x 7 



2 I Q^em 



2tt 



in 



+ £n(x 2 (l - x 7 )) 



^ (1 - * 7 )* 2 /4 J 



+ (l - 3 cos 2 0. 



2(1 -x- 



(86) 



If we integrate Eq. ( |56"D over 7 , we obtain the analogous result derived previously in Ref. [|J. 

For the one-loop fragmentation contributions from parton c(= g, q, q) to 7, the subtrac- 
tion terms again vanish if x c > x max (z, 5, eh), defined in Eq. (|62"D; thus, 

da (1)mcl 



, ~ (l)iso 



E„ 



if x c > x max (z, 5, e h ) 



(87) 



<i 3 p c ' d 3 p c 

In Eq. (|8?j), the corresponding inclusive hard parts are given in Eqs. ([771) and (f?9|). 

If 1/(1 + e/j) < x c < x max (z,5,6h), the subtraction terms are finite and are given in 
Eqs. fl60|) and (|73|). Combining Eq. ( |60D with the inclusive contribution, we obtain the gluon 
fragmentation contribution as 

, ~ (l)iso 



E n 



e+e~->gX 



d 3 Pg 



2E 



a 2 N - 



X 



9. 



c f 



x <M1 + cos 2 7 ) 



1 + (1 - Xg 



.1, (■ 



2tt 



in 



+ £n 



J^\+in(x 2 g (l-x g )) 

"(1 + e fe )(x 7 /x g ) - l/x g 
(1 - x g )5 2 /4 



+ (1 - 3 cos 2 O 



2(1 - x 5 



(88) 



for 1/(1 + e^) < x 3 < x max (,2, 5, e^). The quark fragmentation contribution is 



da 



(l)jso 
e+e-^qX 



d 3 Pq 



2 A/" _L 



x (1 + cos 2 61, 



1 + x 2 

1 — Xi 



+ in (x\{l - Xi)) 



(1 + e h )(x J /x 1 ) - 1/xi 



(1 - Xl )5 2 /4 



2 Vl-xi 



(3xi — 5) 



1 - 3 cos 2 7 ) I 



(89) 
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for 1/(1 + e h ) <xi< x max (>, 5, e h ); x max (z, 5, e h ) is defined in Eq. (||). 

It is a common feature of Eqs. ([86]), ( |88| ) and fl39|) that all three develop a logarithmic 
divergence ^n(l + e/j — l/x 7 ) as i 7 — > 1/(1 + e^). This logarithm results from the collinear 
singularities of the Feynman diagrams shown in Figs. ^, [7] and |8|a. It corresponds to the 
situation in which an unobserved parton inside the isolation cone becomes almost collinear 
with the observed parton. Since x 7 > 1/(1 + e^), which is the same as x g or x q > 1/(1 + 
emin(z)) for all z, the counter terms, defined through Cg> in Eqs. (|2ll), (p6|) , and (|27|) , vanish. 
Consequently, the collinear singularities have no corresponding canceling counter terms. 

This problem is caused by the incompatibility between collinear fragmentation, which 
was used to define the counter terms in Eqs. ([H]), ( p6|) and (|27D, and the cone fragmentation 
used to define the partonic cross section E 7 da^- b -+ c x/ 'd 3 £ with c = 7,?,?, g in Eqs. fl2T|) 
(p6|), and (^7]). As we pointed out in Section |, to deal with the region x 7 > 1/(1 + e/J 
it may be necessary to revise our concept of photon fragmentation functions in the case of 
isolated photons. 



D. One-Loop Contribution to the Isolated Cross Section when x 7 = 1/(1 + e^j 

As x 7 approaches 1/(1 + e/J, a number of the expressions derived above for the isolated 
cross section develop a logarithmic divergence of the form £n|(l + e^) — l/x 7 |. However, 
the value of is an arbitrary parameter chosen in individual experiments. Certainly, a 
completely consistent theoretical prediction should not be sensitive to an arbitrary exper- 
imental parameter, such as e^. We argue below that the reason for the unstable result as 
x 1 approaches 1/(1 + e/J is the breakdown || of the conventional perturbative factorization 
theorem for the cross section of isolated photons in the neighborhood of x 1 = 1/(1 + e^). 

In this subsection, we first use our results of the last two subsections to explain the origin 
of this logarithmic divergence. Then, we present a derivation of the isolated cross section 
for x 1 near the value 1/(1 + e^), and we show that the key issues are the isolation condition 
and the finite cone size for fragmentation. 
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We start with the situation in which x 7 approaches 1/(1 + e^) from below. When x 7 is 
less than 1/(1 + €h), the expression for direct production, given in Eq. (|30D, is well-behaved 
as a; 7 approaches 1/(1 + e^). Similarly, the gluonic fragmentation contribution, defined by 
Eqs. fl81|) and fl82|), is also well-behaved as x 1 approaches 1/(1 + e/J. Although x g can equal 
1 when x 1 = 1/(1 + e/,), the £n(l — x g ) term in Eq. ( |5^D gives a finite contribution after the 
integration over z. Thus, Eqs. (|80f) and fl3~Tl) should also be valid at x 7 = 1/(1 + e^). 

On the other hand, the quark (or antiquark) fragmentation contribution develops a 
logarithmic diverg approaches 1/(1 + e^). This divergence is caused by the 1/(1 — 

Xi) terms in Eq. fl84D , and it can be understood as follows. Consider the following general 
integral 

I{x 1: e h )= f —{t— — ) ^ m ( x -x 1 )F(z,x 1 = x 7 /z), for x 7 < 1/(1 + e h ), 

" , l/(l+e h ) Z \l — X\J 

(90) 

where m = 0, 1, and F(z,x\ = x y / z) is any smooth function over the region of integra- 
tion. The integral 7(x 7 , e^) can be thought of as a simplified version of the one-loop quark 
fragmentation contribution defined in Eq. ([83]). The factors 1/(1 — X\) and £n(l — x{) are 
typical of the short- distance hard part. It is straightforward to perform the integration, and 
we find 

I(xry,e h ) = / )£n rn (l - xi)F(z = x 1 /x 1 ,xi) 

Jx 1 X\ \ 1 — Xi / 

-tn m+l (1 -x 7 (l + e h )) ^±oo as x 7 -f 1/(1 + e h ). (91) 

This example shows that 1/(1 — X\) terms in the short- distance hard part make the isolated 
cross section very sensitive to the value of e^, as x 7 approaches l/(l + e/ l ). The perturbatively 
calculated cross section for isolated photons becomes logarithmically divergent at a different 
value of x 7 if one chooses a different value of e^. 

The reason for this unsatisfactory logarithmic sensitivity is the infrared divergence asso- 
ciated with the limit in which the final-state gluon's momentum goes to zero, which is the 
same as x\ — > 1. Normally, such an infrared divergence is canceled by virtual diagrams. For 
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example, in the calculation of the inclusive contribution, the infrared singularity associated 
with Xi — > 1 from the real gluon emission diagrams, sketched in Fig. |8|a, is canceled by 
the infrared contribution from the virtual gluon exchange diagrams, sketched in Fig. ^b, 
proportional to 5(1 — xi). However, as we demonstrate below, perfect cancellation between 
real and virtual diagrams is upset by the isolation conditions. 

Isolated photons are defined to be photons accompanied by less than E max = etE^ of 
hadronic energy in the isolation cone. For photons produced from parton fragmentation, 
as sketched in Fig. a situation can arise in which E max in the isolation cone is provided 
completely by parton fragmentation. Consequently, no other soft gluons are allowed to enter 
the isolation cone. The phase space for the final-state non-fragmenting partons is reduced 
by the size of the cone. Due to the mismatch of the phase space, the infrared divergences 
associated with the final-state real gluons cannot be completely canceled by the virtual 
diagrams. Such uncanceled infrared contributions vanish if the cone size is zero. 

As an example, we use the one-loop quark fragmentation contribution to the cross section 
for isolated photons to demonstrate the breakdown of the perfect cancellation of infrared 
divergences. In order to make our presentation parallel the statements of the last paragraph, 
we calculate the isolated cross section from the quark fragmentation term directly without 
going through the subtraction term ||. 

To examine the situation in which x 7 approaches 1/(1 + e/J, we must evaluate both the 
real and the virtual diagrams shown in Fig. || The squared matrix element for the real 
diagrams is obtained from Eq. (|65|): 



real 



(e//) 2 (#//) 2 j (l + cos 2 ^ -2e 
+ (l + cos 2 6i - 2e 
+ (l - 3 cos 2 6 1 



Xi 



X\ \ X\ 



2/13 



1 - Xl J 2/13 l-xi 



2/13 



+ 



2/13 



+ 2 



X\ 



(92) 



Equation (|92|) holds for both the inclusive and the isolated cross sections. The key difference 
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resides in the limits of integration over y± 3 = yi 3 /xi. Integrating $13 from to 1, we obtain 
the real contribution to the inclusive cross section; all infrared divergences associated with 
X\ — > 1 are canceled by the contribution from the virtual diagrams On the other hand, 
the isolation conditions require that the 7/13 integration be divided into three regions: 

rvain[y c ,y m ] rl-y~c rl 

#13 =>• / #13 + / dy 13 + / dy 13 , (93) 



10 Jy c Jmax[(l-i/ c ),(l-jM»)] 

where y c and y m are defined in Eq. fl67l). In the first region, the condition < yi 3 < y c ensures 
that a gluon is in the isolation cone of the fragmenting quark; and condition y 13 < min[?/ c , y m ] 
ensures that the total hadronic energy in the isolation cone is less than E max = ey l E 1 . 
Similarly, the condition max[(l — y c ), (1 — y m )} < y 13 < 1 for the third region ensures that 
the antiquark is in the isolation cone and that the total hadronic energy in the isolation cone 
is less than E max . The second interval represents the situation in which neither the gluon 
nor the antiquark is in the isolation cone. 

If y c < l/m, Eq- ( p3| ) shows that the isolated cross section is the same as the inclusive cross 
section, consistent with Eq. (pHp. When y c > y m , (i.e., (1 — £i)<5 2 /4 > z[(l + e^) — l/x 7 ]), 
Eq. fl9"3] ) states that the phase space of the final state gluon (and/or antiquark) is smaller 
than in the case of the inclusive cross section. 

If x-y < 1/(1 + qJ, or, equivalently, y m = z[(l + eh) — l/^ 7 ] < 0, the y\ 3 integration 
defined in Eq. ( p3|) reduces to the second region only 

f l-(l-X!)S 2 /4 

#13^ / #13, (94) 

J(l-si)<5 2 /4 

where we expand y c to order 5 2 . We rewrite the limits of the integration over yi 3 in Eq. (|94"D 

as 

yl-(l-a;i)5 2 /4 rl y(l-a;i)5 2 /4 rl 

/ #13 = / #13 - / #13 - / #13- (95) 

J(l-a:i)<5 2 /4 JO JO J l-{l- Xl )S 2 /A 

Combining the three-particle phase space of Eq. (Q) and the squared matrix element of 
Eq. (^), we derive that the first term on the right-hand-side of Eq. fl9"5D provides the 
complete real gluon emission contribution for the inclusive cross section fl[] 
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(96) 



The superscript (Rx) stands for the contribution of real gluon emission from the first term 
on the right-hand-side of Eq. . In deriving Eq. (PJBJ) , we use 
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The virtual gluon exchange contribution from the diagrams in Fig. |8|b is 

47T/1 2 V 1 1 



(97) 



d 3 px 



at N 
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2^ 



'47T// 2 
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2 3 

7T' 



(98) 

where the superscript (V) stands for the virtual gluon exchange contribution. 

Combining the virtual (V) contribution and a part of the real contribution (Rx), we 
derive 
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+ 
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£n(l — X\) 



1 — X\ 
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2 VI 



X\. 



+5(1 - Xl ) f^p-f) -i(3a:i-5) 



(l -3cos 2 #i) | 



(99) 



As expected, Eq. ( |9S|) shows that, except for the 1/e term associated with the collinear 
singularity between the fragmenting quark and the real gluon, the contribution from the 
first term on the right-hand-side of Eq. ( p5[ ) cancels all infrared singularities from the virtual 
diagrams in Fig. |]b. 

The third term in Eq. fl9~5|) does not provide any terms singular in 1/e; it yields some 
finite terms that vanish as 5 2 — » 0, 



da {R3) 



(100) 



However, the second term generates a number of terms with 1/e poles. Neglecting all terms 
of 0(5 2 ) or higher, we obtain 
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(101) 
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Collecting the contributions in Eqs. (|99"D, ( |100| ), and (|101|) , we find the complete next-to- 
leading order isolated partonic cross section for e + + e~ — > q{pi) + X, 



Ey 



da 



(l)iso 
e+e~— >qX 
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(102) 



As indicated in Eq. (|94"D, i/13 is always larger than zero for a fixed value of xi ^ 1. Therefore, 
Eida^+™° x /d z py defined in Eq. ( |102| ), is independent of ,a|jg. Furthermore, because there 
are no collinear subtraction terms for the isolated partonic cross section, the short-distance 
isolated partonic cross section a^ lso = a^ tso given in Eq. ( |102| ). We note that for xi ^ 1, 
Eq. (|102|) reduces to Eq. (0), as it should. However, as Xi — > 1, the uncanceled 1/e 2 
and 1/e poles in Eq. (|102|) for gnal a breakdown of conventional perturbative QCD 

factorization ||. 

The singularities in Eq. (|102|) corresponding to the uncanceled poles are infrared in 
nature and, as expected, are proportional to 5(1 — xi). As explained above, the uncanceled 
poles come from the interval specified by the second term in Eq. fl95|). This interval results 
from the restricted phase space for soft real gluon emission due to the isolation constraints. 
These poles would be irrelevant if xi ^ 1. However, xi = x^/z, and xi = 1 is kinematically 
allowed when x 7 = 1/(1 + e^). 

We conclude that the conventional factorization theorem for the cross section of isolated 
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photons in e + e annihilation necessarily breaks down when x 7 ~ 1/(1 + e^). 

V. NUMERICAL RESULTS AND DISCUSSION 

In this section, we present numerical results for the isolated photon cross section in 
hadronic final states of e + e~ annihilation, and we discuss the consequences of the breakdown 
of conventional perturbative factorization. 



In the analysis of Section [TV], three regions are identified. In the first region, a; 7 < 
1/(1 + e/j), the isolated cross section is well behaved in perturbation theory except for an 
infrared logarithmic divergence as x 1 approaches 1/(1 + e^), associated uniquely with the 
0(a s ) quark-to-photon fragmentation component of the cross section. The second region 
is the singular point x 7 = 1/(1 + e/J. At this point, the isolated cross section is undefined 
in perturbation theory. It is the point at which the conventional factorization theorem 
breaks down. The third region is the region x 7 > 1/(1 + e^). In this third region, but 
for the incompatibility of the collinear and cone definition of fragmentation, there should 
be no subtraction term and the isolated and inclusive cross sections would be identical. 
However, because of this incompatibility, the 0(a s ) subtraction terms are finite in a very 
small interval 1/(1 + e^) < x 7 < x™ ax , and the perturbative cross section is predicted to 
show a logarithmic divergence as x 7 approaches 1/(1 + e^) from above. This logarithmic 
divergence has a collinear origin. For x 7 > x™ ax the isolated cross section equals the full 
inclusive cross section. All of these features are illustrated in the numerical results presented 
below. 

The analytical expressions, Eqs. fl30|), (|81|) and (|33"D, for the cross sections with x 7 less 
than 1/(1 + eh) should be useful for comparison with data from LEP. The expressions for 
x™ ax > x 7 > 1/(1 + Eqs. (|86|), (|S8|) and ( p9|) should not be taken seriously because of the 
logarithm, £n((l + e^) — l/a; 7 ), caused by the incompatibility between the two definitions of 
the parton fragmentation. For x 7 > x™ ax , the isolated cross section equals to the inclusive 
cross section. Although the phase space above x 7 = 1/(1 + e/J is small, the behavior of 
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the expressions for x 7 > 1/(1 + e/J has important implications for the isolated photon cross 
section in hadronic collisions. In hadronic collisions, the influence of this limited phase 
space is enhanced after the hard cross section is convoluted with beam hadrons' parton 
distributions. 



A. Numerical Results 

Using the analytical results derived in the previous section, we evaluate the cross section 
for isolated photons at the LEP energy yfs ~ Mz, and we compare this cross section 
with corresponding cross section for inclusive photons. In deriving our numerical results, 
we set Mz = 91.187 GeV and Tz = 2.491 GeV. The vector and axial- vector couplings 
and other constants are taken from Ref. The weak mixing angle sin 2 ^ = 0.2319. 



For the electromagnetic coupling strength a em , we use the solution of the first order QED 
renormalization group equation, and let a em (M z ) = 1/128 For the 0(a s ) contributions, 
we employ a two-loop expression for a s (fi 2 ) with quark threshold effects handled properly. 
We set A ( q ] cd = 0.231 GeV, which corresponds to a s (M z ) = 0.112. 

For the quark-to-photon and gluon-to-photon fragmentation functions, we use the an- 
alytical expressions provided in Ref. [16||. These fragmentation functions are leading-order 



functions. In principle, it would be more consistent to use MS fragmentation functions 
evolved in next-to-leading order |l[ . However, since our primary purpose here is to provide a 
theoretical framework for understanding the behavior of isolated photon production, not nec- 
essarily to present the most up-to-date numerical predictions, we believe these leading-order 
fragmentation functions are adequate for demonstrating the features of the cross section for 
isolated photons. 

In presenting results, we divide our cross sections by an energy dependent cross section 
o"o that specifies the leading-order total hadronic event rate at each value of \fs: 

"2 
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This procedure divides out overall multiplicative factors, and, in doing so, shows approxi- 
mately what fraction of the total hadronic rate is represented by prompt photon production. 
One could, of course, divide instead by the next-to-leading-order total hadronic event rate, 
differing from Eq. ( |103| ) only by the overall factor (1 + cis/n). 

In Fig. [TTJ, we compare the cross sections for isolated and inclusive photons as a function 
of the photon energy E 1 at yfs =91 GeV and # 7 = 90°. We set the renormalization scale 
equal to the fragmentation scale, and both of them equal to photon's energy. The isolation 
cone size is chosen to be 5 = 20°. The isolation energy parameters are chosen as = 0.05 
and 6h = 0.15 in Fig. |ll|(a) and Fig. [Tl|(b) , respectively. The value = 0.05 is typical of 
LEP experiments, while = 0.15 is used by CDF and DO at the Fermilab Tevatron collider. 
As expected, the isolated cross section is finite and well-behaved in most of interval of E 1 
less than the critical value \fs /2(1 + e h ) (or i 7 <1/(1 + e^). The perturbatively calculated 
isolated cross section is less than the inclusive cross section in most of this interval, as it ought 
to be. When E 1 approaches y / s/2(H-e /l ) (or a; 7 — > 1/(1 + ^)), the perturbatively calculated 
isolated cross section becomes larger than inclusive cross section, which does not make 
physical sense, and it eventually approaches infinity when x 7 = 1/(1+6^), where conventional 
factorization breaks down. When E 1 approaches the critical value from above (i.e., x 1 > 
1/(1 + 6h)), the perturbatively calculated isolated cross section approaches negative infinity, 
a result of the term proportional to £n((l + e^) — l/a; 7 ) in Eqs. (|86|) , (|88|) and (|89|). 

Figure |Tl~|(b) shows that the perturbatively calculated isolated cross section starts at 
negative infinity when E n = y / i/2(l + e/J. It increases as E 1 increases and quickly be- 
comes equal to the inclusive cross section. Since the values of x^ ax (5, e/J in Eq. (|50| ) and 
x max (z,5,eh) in Eq. fl62|) are very close to the critical value 1/(1 + e/J, the region in which 
Eqs. ([86]), (|88| ) and ( |89| ) are applicable is extremely small. Therefore, the isolated cross 
section is equal to the inclusive cross section in most of the interval of E 1 larger than the 
critical value a/s/2(1 + e^) (or x 1 > 1/(1 + e^)). 

In Fig. [12], we show the breakdown of the isolated cross section in terms of contribu- 
tions from individual pieces. The isolated cross section is plotted as a function of E 1 at 
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y/s = 91 GeV and scattering angle 9 1 = 90°. We choose 5 = 20°, e h = 0.15 and e h = 0.30 
in Fig. |12"Ka) and Fig. |T2](b), respectively. We use the label "0th- Frag" for the lowest-order 
fragmentation contribution, "Direct" for 0(a em ) direct production, and "q-Frag" and "g- 
Frag" for the 0(a s ) quark and gluon fragmentation contributions, respectively. The curves 



in Fig. show that the breakdown of factorization is associated with the 0(a s ) quark frag- 
mentation contribution, which is infrared in nature. All one-loop contributions approach 
negative infinity as E 1 approaches the critical value from above, confirming that this di- 
vergence is of collinear origin, caused by the incompatibility between the two definitions of 
parton fragmentation: exact collinear long-distance fragmentation and the short- distance 
contribution with a finite cone size. 

In Fig. [13|, we examine the predicted # 7 dependence of the isolated cross section when 
Ky is about half of the critical value. We plot the isolated cross section as a function 
of scattering angle # 7 for E 7 = 20 GeV. Other parameters are chosen as shown in the 
figure. There is no contribution from the leading order quark fragmentation, and the 0(a s ) 
contributions are negligible at this energy. The angular dependence of the isolated cross 
section at y/s ~ 20 GeV is determined by the direct contribution, which has both 1 + cos 2 # 7 
and sin 2 6 1 contributions (c.f., Eq. 



In Figs. |TTHT3|, we set the renormalization and fragmentation scale to be // = E~. De- 



pendence of the isolated cross section on fi is examined in Fig. [14]. We plot the isolated 
cross section as a function of [i/E^ at y/s = 91 GeV, 9 1 = 90°, and a fixed E^ = 20 GeV. 
Other parameters are as shown in the figure. As expected, the isolated cross section has a 
very small scale dependence. Isolation eliminates most of the fragmentation contributions, 
and the dominant direct contribution, Eq. fl80"D , is not affected by strong interactions. It has 
neither renormalization nor fragmentation scale dependence. 

Dependence on is strong in the region near a; 7 ~ 1/(1 + e^), but the perturbative 



calculation is not reliable in this region. In Fig. 15, we plot the isolated cross section as a 



function of at a fixed value of £L = 20 GeV, about half of the critical value. Relatively 



small dependence on is observed. As expected, Fig. [L5] shows that the isolated cross 
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section grows as increases. At larger e^, more hadronic energy is allowed in the isolation 
cone, and consequently, more events contribute to the isolated cross section. 

The isolated cross section is predicted to show a specific £n5 2 dependence on the cone 
size 5. This prediction is valid for small values of 5 since we neglect throughout terms 
proportional to finite powers of 5. In doing our analysis, we explicitly constrain 5^0. If 
5 is set to zero before the phase space integrals are done (e.g., in Eqs. (|53|), (|69|), and (|95|)) 
then we would find a^ sub = and a tso[ = a incl . For x 7 < 1/(1 + e^), the £nd 2 dependence 
of the isolated cross section is displayed explicitly in Eqs. (|80"D , Cj82|) , and (jS4"t). 



B. Discussion 



According to the conventional factorization theorem of perturbative QCD [|1(|, the one- 



loop fragmentation contributions derived here become the collinear counter-terms for the 
two- loop direct contribution to e + e~ —>■ 7A at order 0(a em a s ). In similar fashion to the 
derivation of Eq. fl2T|), we apply Eq. (Rf) perturbatively to the order a em a s . We obtain 



(2)sub 



E q (orEq)>e h E~, 

+ tTUxi^Dfl^z) + a^t qX (^)®D(%(z) 

+ (q^q) . (104) 

Since the zeroth-order subtraction term ^+5-^ ^ vanishes, and the zeroth order photon- 
to-photon fragmentation function D^ ry (z) — 6(1 — z), the factorized expression for the 
two-loop short- distance hard-part becomes 

*(2)«i6 / x _ (2)sub 
a e+e-^X\ X l) — a e+e-->-yX 



E q (orE q )>e h E^ 

^l q x^)®D^(z) - a^t qX ^)®D^(z) 
~ -T:-t gX (^D^(z) -{q^ q) . (105) 

It is clear from Eq. ( |105| ) that the one-loop subtraction term, &^- b x derived in this paper, 
Eqs. (|71~D and (0), becomes a collinear counter-term for the two-loop direct contribution. 
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Therefore, due to the breakdown of factorization at x 1 = 1/(1 + e^) at one-loop order 
and the fact that the perturbatively calculated isolated cross section is larger than the 
corresponding inclusive cross section, the higher-order contributions to the isolated cross 
sections may be ill-defined as well. To resolve this issue, it will be important and necessary 
to investigate possible alternative definitions of the cross section of isolated photons to all 
orders in perturbation theory. 

In Eq. (Q), we define the isolated cross section as a difference between the inclusive cross 
section and a subtraction term. The inclusive cross section is well-defined in perturbative 
QCD Q. Since the isolated cross section is an experimental measurable and known to be 
finite, the subtraction term, defined in Eq. (0), should be finite as well. It is not difficult to 
show, as we have in this paper, that the subtraction term is an infrared sensitive quantity in 
perturbation theory. Therefore, the issue is whether we can find a scheme to consistently fac- 
torize the infrared sensitive quantity into a perturbatively calculable infrared safe hard-part 
convoluted with a well-defined long-distance universal function, such as the fragmentation 
function. 

Although the breakdown of factorization is demonstrated in this paper only at one-loop 
level in perturbation theory, the following intuitive argument generalizes our thinking. As 
remarked in Sec. |T|, the subtraction term can be viewed as a "cross section" for a photon "jet" 
with photon momentum i and hadronic energy E^ ne in the "jet" cone, with the hadronic 
energy required to be larger than -E max = e^-E^. This "jet cross section" is an integrated jet 
cross section, where "integrated" denotes the cross section for "jet" events with hadronic 
energy greater than a fixed minimum value (-E max — £h,E^) and up to a whatever maximum 
value is allowed by the kinematics. Experience with the definition and calculation of jet cross 
section fll4"|,|rT| suggests that the integrated "jet cross section" could be finite in perturbative 



QCD, as long as we do not specify details within the jet. However, if we insist on determining 
the energy of one specific parton within the "jet" , and ask how this parton fragments into 
the observed photon, we will most likely encounter an incomplete cancellation of infrared 
and collinear singularities. Introduction of parton-to-photon fragmentation functions for the 
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cross section of isolated photons is similar to asking for the details within the photon "jet". 
Further analysis of these issues is beyond the scope of the present paper. 

C. Phenomenology 

The 0(a em ) direct contribution to the isolated cross section is dominant over most of the 
region x 1 < 1/(1 + th) where perturbation theory is reliable. This statement is illustrated in 
Figs. [3~2] - p~5| . Correspondingly, a direct comparison of our expressions with data on isolated 
prompt photon production is tantamount primarily to a verification of the 0(a em ) expression 
presented in Eq. (|80|) . In this expression, we display the explicit dependence of the cross 
section on E 7 , 9 7 , and 5. On the other hand, an important goal of investigations of e + e~ — > 
7X is also the extraction of quark-to-photon and gluon-to-photon fragmentation functions 
fTT|^T3| . Our analysis of this paper indicates that the task is not straightforward. Not only is 
the direct contribution dominant over a large interval if e/j is small, but as we demonstrated 
to order 0(a s ), the perturbatively computed quark-to-parton fragmentation contribution to 
the isolated cross section is infrared singular in the neighborhood of a; 7 = 1/(1 + e^). This 
means that we cannot offer a reliable theoretical expression for comparison with experiment 
in this region. We expect that further theoretical work in the near future will lead to a 
resolution of this unsatisfactory situation. Meanwhile, we suggest a possible procedure for 
analysis of the data. We suggest that the 0(a em ) direct contribution, specified in Eq. (jBOl). 
be subtracted from the data. Since the direct contribution is of purely electromagnetic origin 
at 0(a em ), it is not subject to hadronic uncertainties. After the subtraction is done, one is 
left with a sample of events that may be attributed to 0(a s ) quark-to-photon and gluon-to- 
photon fragmentation. It would be instructive to examine the E 11 7 , and S dependence of 
this subtracted sample. 

Isolation enhances the 0(a em ) direct contribution and eliminates the leading-order quark- 
to-photon fragmentation contribution for x 7 less than the critical value 1/(1 + e^). It has 
been suggested |2j that the region in which the contribution from leading-order quark-to- 
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photon fragmentation is important can be enlarged if is chosen to be relatively large. 
In Figs. 0(b) and |l2|(c), the dashed curve shows the contribution from the leading-order 
fragmentation contribution; it contributes only in the region above x 1 = 1/(1 + eh). The 
leading-order fragmentation contribution dominates in this region, and the isolated cross 
section is essentially the same as the inclusive cross section [I]. Experimental measurement 
of prompt photons in this region is difficult both because the rate is low and because clean 
identification of a prompt photon signal is problematic if €h is relatively large. When is 
large, the hadronic energy in the cone makes it difficult on an event-by-event basis to be 
certain the photon does not result from tc° or rj decay, or from other hadronic processes. 

We turn next to implications of our analysis for the computation of isolated prompt 
photon cross sections at hadron colliders. In hadron reactions, the lowest-order subprocesses 
are the two 0(a s a em ) direct subprocesses, q + g — > q + 7 and q + q — > g + 7, along with the 
lowest-order fragmentation subprocesses. In the lowest-order fragmentation, various 0{a 2 s ) 
two-parton to two-parton subprocesses are followed by leading-order long distance quark- 
to-photon or gluon-to-photon fragmentation. In computing prompt photon production to 
next-to-leading order in QCD, one must compute all 0{a^a em ) direct subprocesses and 
include fragmentation at next-to-leading order. It is at this level that problems of the 
type we find in our study of electron-positron annihilation will also be encountered in the 
hadronic case. Specifically, the next-to-leading order quark fragmentation contribution is 
ill defined owing to the infrared problem identified in our analysis. In their treatment of 
prompt photon production, Baer, Ohnemus, and Owens JF] include the direct contributions 
through 0{a 2 s a. em ), but they include fragmentation through lowest-order only. Their results 
are therefore not affected by the infrared problem. The studies of Gordon, Vogelsang, and 
collaborators |9j include fragmentation through next-to-leading order. 

It is not our intent in this paper to provide a new analysis of the hadronic case. Rather, 
we wish to point out that the infrared sensitivity of the next-to-leading order quark-to- 
photon fragmentation contribution is enhanced in the hadronic case. We also suggest a 
temporary phenomenological remedy. In hadron-hadron collisions, we are interested in the 
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production of prompt photons as a function of the photon's transverse momentum, p?. 
Invoking factorization, and continuing to use I to denote the four-vector momentum of the 
photon, we write the inclusive or the isolated cross section in hadron(A) + hadron(B) — ► 7X 
as 



The sum is taken over all initial partons (q, q, g) in the incident hadrons, and <j>(x) are parton 
density distributions. The partonic level cross section is expressed as 



where c = 7, q, q, or g. If the factorization theorem holds for isolated photons, the partonic 
hard-parts a a b^cx m Eq. ( |107|) should be finite for isolated photons and should be smaller 
than the corresponding hard parts for the inclusive cross section. This is not the case for the 
perturbatively computed quark fragmentation contribution at one-loop, as we now illustrate. 
To streamline the discussion we begin by defining a parton-parton flux $(r), 



where r = x a x^. To exploit our electron-positron results most conveniently, we limit our 
discussion to the contribution in Eq. ( |106|) associated with the quark-antiquark initial state, 
and, further, we take only the term in which the flavors of the initial and final quarks differ: 
q' + q' — > q + q + g, where q fragments to 7. We point out, however, that there is a next- 
to-leading order quark fragmentation contribution associated with almost all subprocesses. 
Our illustration is therefore much more general than the specific example we treat. We 
further specialize to prompt photon production at zero rapidity (i.e., at 90° in the hadron- 
hadron center-of-mass frame), and ignore real gluon emission from the initial-state quark (q') 
antiquark (<f) and from the s-channel virtual gluon (g). As an expedient approximation we 
take equal values for the partonic fractions, x a = Xb = x = yfr. This latter approximation 
allows us, in the isolated case, to use our Eqs. ( j83|) and (|84] ) directly, setting 6 1 = 90° in 




(106) 




(107) 




(108) 
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Eq. fl34|). Likewise, in the inclusive case, we use our Eqs. (|78| ) and (|79|) directly, setting 
6-y = 90° in Eq. (|79|). The normalization factor differs, of course, in the hadron case. In 
Eqs. (|79| ) and (|84"D, we replace 

s — > s = x a Xf,s = t s ; (I09a) 

-< c (*) - 4 ; ( 109b ) 

«L ~> " s 2 ; (109c) 
iV 2 - I 

N e C F ^^-. (109d) 



In the translation to the hadronic case, the variable x 7 becomes xt where xt = 2pt/\s ~ 
xt/x with xt = 2p T / v ^s. 

The special one-loop quark fragmentation contribution to the observed cross section takes 
the form 

E ^ aA ^ X ~ dr <$> q , q ,{r) E^ d(Tq ' q '^ X (r) + other subprocesses . (110) 



In Eq. ( |110| ), the one- loop contribution E^da q >qi-+ 7 x/d 3 £( T ) ^° inclusive [or isolated] cross 



section is obtained from Eq. ( |75|) [or Eq. (|54])1 with the replacements defined in Eq. 



In Fig. [16|, we present the one-loop quark-to-photon fragmentation contributions to the 
inclusive and isolated prompt photon cross sections as a function of the scaling variable x 1 in 
e + e~ annihilation. The region of infrared instability in the neighborhood of x 7 = 1/(1 + e^) 
is evident, as is the fact that the isolated cross section exceeds the inclusive in part of the 
Xry interval. 

In Fig. [17|, we show the integrand in Eq. (|110|) obtained after convoluting the parton flux 



of Eq. (|108|) with the one- loop contributions in Fig. [16] (with the replacements in Eq. ( |109| )) 



In obtaining Fig. [L7], we use CTEQ3M parton distributions |18[ to compute the effective 



parton flux, $ g /g/(r) with q' = u (up quark); the final-state quark flavor is chosen as q = d 
(down quark). The features demonstrated in Fig. [17] are independent of the choice of quark 
flavors. 
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The one-loop fragmentation contribution to the hadronic cross section at a given value of 
transverse momentum is obtained by computing the area under the curve in Fig. [17] from 
xt = 2px/\fs to 1. It is evident that the convolution with the parton flux substantially 
enhances the influence of the region of infrared sensitivity. The divergences above and 
below the point xt = xt/\/t = 1/(1 + e^) [or x = y/r = xt(1 + e h)] are integrable 
logarithmic divergences, and thus they yield a finite contribution if an integral is done over 
all x (or all r). However, we stress that the perturbatively calculated one- loop partonic cross 
section E q da q 1?,^ qX , obtained from Eq. (]84D , has an inverse-power divergence as x q — > 1 
and has uncanceled 1/e poles in dimensional regularization ||. The divergence for a; 7 < 
1/(1 + €h) becomes a logarithmic divergence only after convolution with a long-distance 
fragmentation function. The theory is not defined at the point Xt = 1/(1 + 6^), and it would 
be inappropriate to take the perturbatively calculated cross section at face value when it 
exceeds the inclusive cross section. 

As remarked earlier, a new theoretical definition of the isolated photon cross section is 
desirable, along the lines of a "photon jet" definition. It would be infrared safe to all orders in 
a s . Presumably, the new definition would avoid parton-to-photon fragmentation functions. 
Meanwhile, we propose a simple phenomenological remedy to deal with infrared problem 



illustrated in Figs. [16] and [17J. In calculations of isolated prompt photon cross sections in 
hadron-hadron reactions, we suggest that the perturbatively calculated one-loop quark-to- 
photon fragmentation contribution to the isolated cross section be used only in the region 
of phase space in which it is calculated to be smaller than the inclusive rate. Elsewhere, the 
one-loop inclusive contribution should be used as an upper bound. 
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FIGURES 



FIG. 1. Illustration of an isolation cone about the direction of the photon's momentum. 
FIG. 2. Illustration of e + e~ — > jX in an m parton state: e + e~ — ► followed by fragmentation 

FIG. 3. Illustration of an isolation cone containing a parton c that fragments into a 7 plus 
hadronic energy Ef rag . In addition, the cone includes a gluon that fragments giving hadronic 
energy E parton . 

FIG. 4. Lowest order, O (a° m a°), photon production through quark fragmentation. 
FIG. 5. Feynman diagrams for e + e~ — > 799. 

FIG. 6. Order a s Feynman diagram for the parton level fragmentation function D q 1 X-y{z). 

FIG. 7. Order a s Feynman diagrams for e + e~ — > 555 that contribute to e + e~ — ► 7 X via 5^7 
fragmentation. 

FIG. 8. Contributions to the 0(ct s ) cross section ^g+g-^x' both real gluon emission diagrams 
(e + e _ — > and virtual gluon exchange terms are drawn. 

FIG. 9. Order a s Feynman diagrams for the parton level fragmentation function Dq~\ q (z). 

FIG. 10. Center of mass coordinate axes of an e + e~ collision with the z-axis being the direction 
of the observed photon. 

FIG. 11. Comparison of the normalized cross sections for isolated and inclusive photons in 
e + e~ — > 7X as a function of photon energy Ey. (a) = 0.05; (b) = 0.15. 
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FIG. 12. Normalized isolated photon cross section as a function of photon energy E^ : 
(a) = 0.15; (b) = 0.30. The solid line is a the sum of contributions from all individ- 
ual subprocesses. Shown also are the individual contributions from the separate subprocesses. 
The 0(a em ) direct contribution (upper dotted line) nearly saturates the total contribution for 
x 7 < 1/(1 + eh). The dashed line for x > 1/(1 + e^) shows the 0(aP em a° s ) fragmentation contribu- 
tion. The 0{a s ) quark-to-photon and gluon-to-photon fragmentation contributions are illustrated 
with the dot-dashed and the lower dotted lines, respectively. 

FIG. 13. Normalized isolated photon cross section at E 1 = 20 GeV as a function of 9 1 . The 
solid line is the sum of all contributions. The 0(a em ) direct contribution (dashed line) essen- 
tially saturates the total. The 0(a s ) quark-to-photon (dot-dashed) and gluon-to-photon (dotted) 
contributions are also shown. 

FIG. 14. Normalized isolated photon cross section as a function of n/E y . There is no leading 
order fragmentation contribution: "0-Frag" at E 1 = 20 GeV. The labeling of the curves is the 
same as in Fig. 13. 

FIG. 15. Normalized isolated photon cross section as a function of e^. There is no leading 
order fragmentation contribution: "0-Frag" at = 20 GeV. The labeling of the curves is the 
same as in Fig. 13. 

FIG. 16. Comparison of the one-loop quark fragmentation contributions to the isolated cross 
section and the inclusive cross section in e + e~ — > jX as a function of x 7 = 2£ l 7 /- v /i. 

FIG. 17. Integrand of Eq. (110) as a function of xt/^/t = xt/x for pp — > at yfs = 1.8 TeV 
and pTp = 20 GeV and rapidity zero. The isolation parameters are = 0.15 and 5 = 20°. 



54 




Fig.l 




Fig.2 



7 




X 
Fig.3 




Fig.4 



7 



Fig 




Fig.6 



9 



Fig 



+ 



(a) 



+ 



(b) 
Fig.8 




Fig.10 



Photon Cross Sections 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

\ 1 1 1 1 


i i i 


~ \ 


- 


r V 






\ \ 






- 


\ x Inclusive 




>v \ 




Isolated \. s \ 

\. \ 

N. v. 


- 


— \ 

! 6=20° x 




e h = 0.05 


\ : 


Vs~ = 91 GeV 




6 7 -90 Q 




- M=E 7 
i i i i 1 i i i i 1 i i i i 1 i i i i 1 i 


i i i 



10 20 30 40 50 



E ? (aeV) 



Fig.ll-(a) 



Photon Cross Sections 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

- v 1 

~ \ 
\ \ 

- V 

- \ x 

\ x 
\ \ 

\ x Inclusive 


1 i i i i 

- 


- 


\ \ 

\ N 

Isolated \. \ 

\. \ 


- 


• (5-20° 


\ 


: e h = 0.15 




Vs~ = 91 GeV 




6 7 -90 Q 




- M=E 7 
i i i i 1 i i i i 1 i i i i 1 i i i i 


1 i i i i 



10 20 30 40 50 



E ? (aeV) 



Fig.ll-(b) 



Isolated Cross Sections 



1 


1 1 1 1 1 1 1 1 1 


i i i i 1 i i i i 
Oth-Frag 




— 




Total 


— 


— 




Direct '.' 

/ 

/ 

s 




z 




q-Frag 


z 




(5 = 20° 








e h = 0.15 
d y =90° 
Vs=91 GeV 


g-Frag 




- 1 


i i i 1 i i i i 1 


i i i i 1 i i i i 


1 i i i - 



10 20 30 40 



E ? (GeV) 



Fig.l2-(a) 



Photon Cross Sections 



10 



-5 



cv 

I io- 



CO 

T3 



10 



10" 



o 
b 



10 



7 



8 



9 



10 



-10 



lo- 



ll _ 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

r \ : 0th Frag "1 


r X^Total 




Direct ' < 




= / 








e~ q-Frag 




_ g-Frag 

<5 = 20° 
: c: h = 0.3 .... 
^ 7 = 9O° 

j V^ = 91 GeV 


- 

• Direct _ 


~ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


i i 1 i i i ~ 



10 20 30 40 
E (GeV) 



50 



Fig.l2-(b) 





1 o w 1 CX U C ^ w 1 Uoo k_) C v_. LI U 1 1 o 


lO" 6 


1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 1 1 1 

r Total - 


cv 

> 

O 

C 10 -7 


Direct 


A _ on 
— cU 


e h = 0.15 




: E 7 20GcV : 


CO 


; Vs=9i Gev ; 


b R 


, , TP 






o 
b 


q-Frag 


" lO" 9 


: g-Frag : 


io- 10 


i i i i 1 i i i i 1 i i i i 1 i i i i 1 i i i i 



20 40 60 80 100 
# 7 (Degrees) 
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